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INTRODUCTION 



! Entanglement is one of the most striking features of quantum theory and is of capital 

\ importance for the whole field of quantum information theory (see, e.g., Refs. The 
determination whether a given quantum state is entangled or separable is still an open and 
f — . challenging problem, in particular for higher dimensional systems. 

! For a two-qubit system the geometric structure of the entangled and separable states in 

\ the Hilbert-Schmidt space is very well known. Due to the Peres-Horodecki criterion 0, 0] 
• we know necessary and sufficient conditions for separability. This case is, however, due to 
O ■ its high symmetry quite unique and even misleading for conclusions in higher dimensions. 

In higher dimensions, the geometric structure of the states is much more complicated 
and new phenomena like bound entanglement occur [g], 0, Isl, 0, 10, 11|. Although we do not 



' know necessary and sufficient conditions a la Peres-Horodecki we can construct an operator 
^ - entanglement witness - which provides via an inequality a criterion for the entanglement 

of the state 0, E, IH, [11 , 



In this paper, we use entanglement witnesses to explore the geometric structure of the 
two-qutrit states in Hilbert-Schmidt space, i.e. we geometrically quantify entanglement for 
special cases, and present a method to detect bound entanglement. Two qutrits are states on 
the 3x3 dimensional Hilbert space of bipartite quantum systems. In analogy to the familiar 
two-qubit case, which we discuss at the beginning, we introduce a two- and three-parameter 
family of two-qutrit states which are part of the magic simplex of states [ll|, 15, 1^ and 



determine geometrical properties of the states: For the two-parameter family we quantify the 
entanglement via the Hilbert-Schmidt measure, for the three-parameter family we discover 
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bound entangled states in addition to known ones in the simplex j^, 1^. Finally, we give 
a sketch of how to use our method to construct the shape of the separable states for the 
three-parameter family. 



II. WEYL OPERATOR BASIS 

As standard matrix basis we consider the standard matrices, the d x d matrices, that 
have only one entry 1 and the other entries and form an orthonormal basis of the Hilbert- 
Schmidt space, which is the space of operators that act on the states of the Hilbert space 
T-C^ of dimension d. We write these matrices shortly as operators 

\j){k\, with j,k = l,...,d, (1) 

where the matrix representation can be easily obtained in the standard vector basis {\i)}- 
Any matrix can easily be decomposed into a linear combination of matrices ([1]). 

The Weyl operator basis (WOB) of the Hilbert-Schmidt space of dimension d consists of 
the following d"^ operators (see Ref. [l3]) 

d-l 

Unm = ^e^*="|fc)((A; + m)modd| n,m = 0,1, ... ,d - 1 . (2) 

k=0 



These operators have been frequently used in the literature (see e.g. Ref s. | 111. Il5l . Il8l. Il9|). 
in particular, to create a basis of maximally entangled qudit states ISl, 120, l2l| . 



Example. In the case of qutrits, i.e. of a three-dimensional Hilbert space, the Weyl 
operators ([2]) have the following matrix form: 



[/oo = 1 , [/oi = 1 , f/o2 = I 1 I , (3) 

\0 1/ \i 0/ 

/lO 0\ /010\ 

f/io = e2-V3 \ , Un= \ e^-^/^ , Uu = 




^^20 = e-^-*/^ , f/21 = e-^-*/^ , U22 





Using the WOB we can decompose quite generally any density matrix in form of a vector, 
called Block vector 17 



d-l ^ ^ 

- 1 + ^ bnmUnm = ^ 1 + 6 " t/ , (4) 

n,m=0 



with n,m = 0,1, ...,d — 1 (600 = 0). The components of the Bloch vector h = ({&„m}) 
are ordered and given by bnm = Tr Unm P ■ In general, the components bnm are complex 
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since the Weyl operators are not Hermitian and the complex conjugates fulfil the relation 
^nm = 6 ^"'"^ b-n-m , which foUows casily from definition ([2]) together with the hermiticity 
of p . Note that for d > 3 not any vector b of complex components is a Bloch vector, i.e. a 
quantum state (details can be found in Ref. [l3])- 



The standard matrices ([T]), on the other hand, can be expressed by the WOB as [17 



1=0 

It allows us to decompose any two-qudit state of a bipartite system in the dx d dimensional 
Hilbert space Ti.^ (S> Ti.^ into the WOB. In particular, we consider the isotropic two-qudit 
state pL'^'' which is defined as follows (g], 10, l^l : 

pL'^ = «K>«| + aeR, -^^<a<l. (6) 

The range of a is determined by the positivity of the state. 
The state 

I J.d 



denotes a Bell state, which is maximally entangled, and in terms of the WOB it is expressed 
by (see Refs. QQ) 



with 

d-l 



U := J^^irn^U.im, (/,m)^(0,0). (9) 



l,m=0 



The negative values of the index / have to be considered as mod d . 
Then the isotropic two-qudit state follows from Eq. ([8]): 



pi^) = + (10) 



We call decomposition (ITU]) Bloch vector form of the density matrix. 



III. ENTANGLEMENT WITNESS HILBERT SCHMIDT MEASURE 

A. Definitions 

The Hilhert- Schmidt (HS) measure of entanglement |23| . |24| is defined as the minimal 
HS distance of an entangled state pent to the set of separable states 5", 

-D(pent) := min Hp - pcntll = ||po-Pcnt|| , (11) 



where po denotes the nearest separable state, the minimum of the HS distance. 
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An entanglement witness A is a Hermitian operator that "detects" the entanglement of 
a state pent via inequahties [l, [H, 13, 14, 17 



(Pent,^) = Trpentv4 < 

{p,A) = TTpA > 



(12) 



An entanglement witness is "optimal" , denoted by Aopt , if apart from Eq. (fT2l) there exists 
a separable state po G S" such that 



(po, Aopt) 







(13) 



The operator Aopt defines a tangent plane to the set of separable states S and can be 
constructed in the following way [ij, l25|] : 



A 



Po - Pont - (Po, Po - Pent) 1 



opt 



I Po ~ Pent I 



(14) 



The term ||po — Pent|| is a normalization factor and is included for a convenient scaling of 
the inequalities (|T2l) only. Thus, the calculation of the optimal entanglement witness Aopt 
to a given entangled state pent reduces to the determination of the nearest separable state 
Po . In special cases, we are able to find po by a guess method fisl] , what we demonstrate 
in this article by working with our Lemmas, but in general its detection is quite a difficult 
task. 



B. Two parameter entangled states — qubits 

It is quite illustrative to start with the familiar case of a two-qubit system to gain 
some intuition of the geometry of the states in the HS space. We aim to determine the 
entanglement witness and the HS measure of entanglement for the following two-qubit 
states which are a particular mixture of the Bell states IV'"), |0^), 10^)^ 

Pa,/5 = 1 + «|0+)(0+| + ^ (|V^+)(V^+| + |r)(n) . (15) 

The states (fT5|) are characterized by the two parameters a and (3 and we will refer to the 
states as the two-parameter states. Of course, the positivity requirement pa,i3 > constrains 
the possible values of a and /5, namely 

a < -f3 + l, a > -f3--, a < (3 + 1, (16) 

3 3 

which geometrically corresponds to a triangle, see Fig. [H 

According to Peres j3] and the Horodeckis (11 the separability of the states is determined 
by the positive partial transposition criterion, which says that a separable state has to stay 
positive under partial transposition (PPT). For dimensions 2x2 and 2x3 the criterion is 
necessary and sufficient |5|], thus any PPT state is separable for these dimensions. States 
(IT^ which are positive under partial transposition have the following constraints: 

a > (3-1, a < ^(3 + ^, a>-(3-l, (17) 
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FIG. 1: Illustration of the two-qubit states pa^i3 (|15p and their partial transpositions 



and correspond to the rotated triangle; then the overlap, a rhombus, represents the separable 
states, see Fig. [H 

In the illustration of Fig. [1] the orthogonal lines are indeed orthogonal in HS space. 
Therefore, the coordinate axes for the parameter a and j3 are necessarily non-orthogonal. 
In particular, the a axis has to be orthogonal to the boundary line a = l/3(/? — 1), and the 
f3 axis has to be orthogonal to a = (3 + 1. 

The two-parameter states po,/3 define a plane in the HS space. It is quite illustrative to 
see how this plane is located in the three-dimensional simplex formed by states that are 
mixtures of the the four two-qubit Bell states. The simplex represents a tetrahedron due to 
the positivity conditions of the density matrix [3, 26, 23|. Applying PPT, the tetrahedron is 



rotated producing an intersection - a double pyramid - which corresponds to the separable 
states. An illustration of the described features is given in Fig. [2l 

To calculate the HS measure (fTTj) for the two-parameter qubit state (fTSll we express the 
state in terms of the Pauli matrix basis 

Pa,p = ^ (1 + a ((Ti (g) (Ti - (72 ® + (a - /?) (T3 (g) 0-3) , (18) 

where we have used the well-known Pauli matrix decomposition of the Bell states (see, e.g., 



Ref. [1|) 



In order to determine the HS measure for the entangled two-parameter states p^"J^ we 
have to find the nearest separable states, which is usually the most difficult task to perform 
in this context. In Ref. [251] a Lemma is presented to check if a particular separable state is 
indeed the nearest separable state to a given entangled one. 

Lemma 1. A state p is equal to the nearest separable state po if and only if the operator 

^ ^ P - Pent - {P,P- Pent) 1 ^^g^ 
Hp - Pent II 



is an entanglement witness. 
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FIG. 2: Location of the plane of states p^^ (jlSp in the tetrahedron formed by the Bell states. 

Lemma [T] probes_ if a guess p is indeed correct for the nearest separable state. If this is 
the case, operator C represents the optimal entanglement witness A^^t ( IT4I1 . 

Lemma [1] is used here in the following way. First, for any entangled two-parameter state 
( ITSl) we calculate the separable state that has the nearest Euclidean distance in the geometric 
picture (Fig. [1]) and call this state p. But since the regarded picture does not represent the 
full state space (e.g., states containing terms like a* 1 or 6j 1 (g) a* are not contained on 
the picture), we have to use Lemma [T] to check if the estimated state p is indeed the nearest 
separable state po- 



1. Region I 

Let us consider first the entangled states located in the triangle region that includes the 
Bell state I*/)"'") , i.e. Region I in Fig. [H An entangled state in Region I is characterized by 
points, i.e. by the parameter pair (a,/5), constrained by 

a < 13 + 1, a < -13 + 1, a > -13 + - . (20) 

3 3 

The point in the separable region of Fig. [1] that is nearest (in the Euclidean sense) to the 
point (a,/3) is given by (| + |/?,/?), which corresponds to the state 

I f 1 + (3 , ,1-2/3 \ 
P0 = ^ I 1 H ^ (CTI ® CTl - 0-2 ® (^2) H — 0-3 ® 0-3 1 • (21) 

For the difference of the "nearest Euclidean separable" and the entangled state, we obtain 

= iC-^-A^, (22) 
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where E is defined by 

E := CTi ® (Ti - (T2 ® 0-2 + 0-3 ® CTs • (23) 
Using the norm ||S|| = 2^/3 we gain the HS distance 

l|p^-pSII = f (24) 

To check whether the state p/3 coincides with the nearest separable state po;/3 in the sense 
of the HS measure of entanglement ( ITTl) (which has to take into account the whole set of 
separable states), we have to test - according to Lemma [1]- whether the operator 

II 7. „ont II V / 

is an entanglement witness. Remember that any entanglement witness A that detects the 
entanglement of a state p^^^ has to satisfy the inequalities (fT2l) . 
We calculate 

(P/3, P> - pTp) = Tr pMp/3 - P^J:^) = (« - ^ - ^Z^) (26) 
and use Eqs. f l22l) and fl24|) to determine the operator C for the considered case 

(27) 

Then we find 

{pT:,.C) = -^[^-l-l^) < 0' (28) 

since the entangled states in the considered Region I satisfy the constraint a > + | . 
Thus, the first condition of inequalities (fT2|) is fulfilled. 

Actually, condition (!28|) is just a consistency check for the correct calculation of operator 
C since by construction of C we always have {p^^^,C) = —\\p — p'^'^^W < 0. Thus more 
important is the test of the second condition of inequalities (fT2ll and in order to do it we 
need the following Lemma. 

Lemma 2. For any Hermitian operator C on a Hilbert space of dimension 2x2 that is of 
the form 

C = a ^1 + ^QCTj (g) o-j^ a G IR+, Cj G R (29) 

the expectation value for all separable states is positive, 

(p,C) > VpG5, if \ci\ < 1 Vi. (30) 

Proof. Any separable state p is a convex combination of product states and thus a separable 
two-qubit state can be written as (see Refs. [3, HI]) 

p = EfcPfc i (1 ® 1 + Ei ® 1 + 1 ® + Eij ^i^'j ® f^i) . 

with nf,mj^ G R, In'"! < 1 , Im'^l < 1 , Pfc > 0, Efci'fc = 1 ) (31) 
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where \n^\'^ = ('^i^)^- Performing the trace, we obtain 



(p,C) = TrpC = E^^^M^ + E' 

k \ i 



,kk 



(32) 



and using the restriction |cj| < 1 \/i we have 



E 



k k 



(33) 



and since the convex sum of positive terms stays positive we get {p,C) > Vp G 5 . □ 

Since the operator C fl27|) is of the form fl29|) we can use Lemma [2] to verify 

(p,C') > WpeS. (34) 

Therefore, C (1271) is indeed an entanglement witness and pp is the nearest separable state 
P/3 = Po;/3 for the entangled states p'^^ in Region I. 



Finally, we find for the HS measure of the states in Region I 

^'(Pa!^) = IIP0;/3 - PT^W 



73 
2 



a B 

3 3^ 



(35) 



Note that the HS measure fl35|) is equal to the maximal violation of the entanglement witness 



inequality (fT^ as it is shown in detail in Refs. [IJ, ll7|, [25 



2. Region II 

It remains to determine the HS measure for the entangled states p^^^j^ located in the 
triangle region that includes the Bell state , i.e. Region H in Fig. [H Here, the entangled 
states are characterized by points (a,/?) , where the parameters are constrained by 



a < p + l, a > -f3--, a < -(3-1. 

3 3 



(36) 



The states in the separable region of Fig. [T]that are nearest to the entangled states (a,/3) 
in Region II are called p^^^ and characterized by the points 



1/3 (-1 + 2a-/3) 
1/3 (-2-2a + /3) 



(37) 



The necessary quantities for calculating the operator C are 



Pa,/3 - PTj^ = (a + 1 + /3) ((Ti ® Ox - 02® O-^® CT3) 



|Pa,/3 - Pa,(i\ 



2^3 



;-a-i-/?) 



(38) 
(39) 
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(P«,/5, Ptp - P-^P) = (« + 1 + /?) , (40) 
so that C is expressed by 

C = (1 + (Ti (g) CTi - (72 ® (72 - (73 (g) (T3) . (41) 

2v 3 

To test C for being an entanglement witness, we need to check the first condition of inequal- 
ities ( IT2l) : we get 

{pl^^,C) = J-(a + l + /5) < (42) 



as expected. Since operator C f HTj) is of the form fl29|) we apply Lemma [2] and obtain for 
the separable states 

{p,C) > VpG5. (43) 

Therefore, also in Region 11, operator C fj4T]) is indeed an entanglement witness and pa,i3 is 
the nearest separable state pa,i3 = Po;ai3 for the entangled states 

For the HS measure of the states in Region 11, we find 

D'\p:^}^) = \\po;a,p - pT}^\\ = ^i-a-l-P) . (44) 

C. Two parameter entangled states — qutrits 

The procedure of determining the geometry of separable and entangled states discussed 
in Sec. 1111 Bl can be generalized to higher dimensions, e.g. for two-qutrit states. Let us 
first notice how to generalize the concept of a maximally entangled Bell basis to higher 
dimensions. A basis of maximally entangled two-qudit states can be attained by starting 
with the maximally entangled qudit state ([7]) and constructing the other (P — 1 states 
in the following way: 

10,) = t = l,2,...,d^-l, (45) 

where {Ut} represents an orthogonal matrix basis of unitary matrices and Uq usually denotes 



the identity matrix 1 (see Refs. |20l. l2l[|). 

A reasonable choice of the basis of unitary matrices is the WOB (see Sec. [Tll)- Such a 
construction has been proposed in Ref. [3]. Then we set up the following (f projectors onto 



the maximally entangled states - the Bell states: 

Pnk := (f/nfc®l)|0t)(0t|(f/lfc®l) n,k = 0,l,...,d-l. (46) 

In case of qutrits {d = 3), mixtures of the nine Bell projectors 0461) form an eight-dimensional 
simplex which is the higher dimensional analogue of the three-dimensional simplex, the 
tetrahedron for qubits, see Fig. [21 This eight-dimensional simplex has a very interesting 



geometry concerning separability and entanglement (see Refs. llSl. Il6l|). Due to its high 
symmetry inside - named therefore the magic simplex by the authors of Ref. 15|] - it is 
enough to consider certain mixtures of Bell states which form equivalent classes concerning 
their geometry. 
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We can express the Bell projectors as Bloch vectors by using the Bloch vector form ([8 
of Pqo •= 10+) (0+1 and the relations (indices have to be taken modd) [l^ 



UL = e^-f/_._, (47) 

UnmUlk = e~"''' Un+lm+k ■ (48) 

It provides for the Bell projectors the Bloch form 

d-l 

d^ 



Pnk = e'5^('='-"'") U,^ ® U.i^ . (49) 



m,l=0 

We are interested in the following two-parameter states of two qutrits as a generalization of 
the qubit case, Eq. (HM . 

Pa,p = ^~ 1 + aPoo + (Pio + P20) . (50) 

According to Ref. [l5| the Bell states represent points in a discrete phase space. The indices 
n, k can be interpreted as "quantized" position coordinate and momentum, respectively. The 
Bell states Pqo, Pio and P20 he on a line in this phase space picture of the maximally entangled 
states, they exhibit the same geometry as other lines since each line can be transformed into 
another one. 

Inserting the Bloch vector form of Pqq, Pio and P20 (HH]) we find the Bloch vector expansion 
of the two-parameter states (|50|) 



Pa,(3 = Ul + (a-^]U^ + {a + P)U2] , (51) 



9 V V 2 
where we defined 

Ui := Uoi ^ Uoi + U02 ^ U02 + Uii ^ U-ii + U12 ^ U-12 + U21 ® U-21 + U22 ^ U-22 , 
U2 := U^ + U^' with := ?7io ® f/-io , U^' := U20 ® U^20 ■ (52) 

The constraints for the positivity requirement {pa,i3 > 0) are 

a <-p+l, a < + «>---, (53) 



and for the PPT 



2' 

a<-/3-l, a>^/?-i, «<^ + l. (54) 

The Euclidean picture representing the HS space geometry of states (150|) is shown in Fig. [31 
The parameter coordinate axes are chosen non-orthogonal such that they become orthogonal 
to the boundary lines of the positivity region, a = f — ^ and a = lf3+l , in order to reproduce 



the symmetry of the magic simplex. It is shown in Ref. [15[ that the PPT states Pa,f3 are all 
separable states, thus there exist no PPT entangled or bound entangled states of the form 
(ISO]) . Bound entanglement we detect for states that need more than two parameters in the 
Bell state expansion (see Sees. IIII D 21 and IIII D"3!) . 

To find the HS measure for the entangled two-parameter two-qutrit states we apply the 
same procedure as in Sec. IIIIBI we determine the states that are the nearest separable ones 
in the Euclidean sense of Fig. [3] and use Lemma [1] to check whether these are indeed the 
nearest separable ones with respect to the whole state space (for other approaches see, e.g.. 



Refs. |28|,l29|) 
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FIG. 3: Illustration of the two-qutrit states Pa,p (|50p and their partial transposition. The regions I 
and II label the regions where the states are entangled, they are PPT and separable in the overlap 
with the region of PPT points. The PPT points become semipositive under partial transposition. 



1. Region I 

First, we consider Region I in Fig. [31 i.e., the triangle region of entangled states around 
the a-axis, constrained by the parameter values 

a < ^-f3+l, a<-f3 + l, «>f + (55) 



In the Euclidean picture, the point that is nearest to point (a,/?) in this region is given by 

+ 



+ 1/3, which corresponds to the separable two-qutrit state 



with Ui and U2 defined in Eq. ( !52|) . 

For the difference of "nearest Euclidean separable" and entangled state, we find 

where U = U1 + U2 (defined in Eq. ([9])). Using for the norm \\U\\ = 3^8 = 6^2 we gain the 
HS distance 

It remains to calculate 

(p/3, Pf3 - Pa,p) = Trpp{pf3 - Pa,f3) = " ^ ~ ^1^^ (59) 



12 



to set up the operator 

g^P>-p5-fe.P.-PS)i^ 1 (60) 

We test now whether it represents an entanglement witness, i.e., whether C fl60l) satisfies 
the inequahties f|T2l) . As expected, we find 

= < 0. (61) 

To check the second condition of inequahties (|T2|) we set up the foUowing Lemma, similar 
to Lemma [21 

Lemma 3. For any Hermitian operator C of a bipartite Hilbert- Schmidt space of dimension 
d X d that is of the form 



C = a I (ci - 1) ld2 + ^ Cnm Unm ® U-nm | , « G IR+, Cnm € C (62) 
\ n,m=0 / 

the expectation value for all separable states is positive, 

{p,C) >OypeS, if \cnm\ < 1 Vn,m. (63) 
Proof. Any bipartite separable state can be decomposed into Weyl operators as 

^ d-l d-l 

P = Xl^fc ^(id ® Id + Yl ^^d^n^^Unm ® Id + XI " 1 ^^Ik Id ® Ulk 

k n,m=0 l,k=0 

di-1 

n,m,l,k=0 

nL,mf, eC, |n^|<l, |m^|<l, Pfc>0, J]pfc = l, (64) 

k 

where we define |n''|^ := J2nm^nmnnm = E„m l^nmP- 

Performing the trace, we obtain (keeping notation formula (l65l) becomes more evident) 



{p,C) = TrpTC = Xpfe|(ci-l)aM +Xc, 
and using the restriction |c„m| < 1 yn,m we have 



(65) 



< Y\<^\\'^-nm\ < 1, (66) 



ra,m 



and since the convex sum of positive terms stays positive we get {p,C) > Vp G 5 . □ 
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Since the operator C (160|) is of the form (162|) we can use Lemma [3] to verify 

{p,C) > ^peS. (67) 

Thus C ( !60l) is indeed an entanglement witness and pp is the nearest separable state 
P/3 = po;/3 for the entangled states p™^ in Region I. 

For the HS measure of the entangled two-parameter two-qutrit states ( l50|) we find 

D\pt,) = Wpo;, - pSII = ^ (" - i - ^/^) • (68) 

Region II 

In Region II of Fig. [3l the entangled two-parameter two-qutrit states are constrained by 
«<^/?-^, «>^/?-^, « < + (69) 
The points that have minimal Euclidean distance to the points (a,/?) located in this region 



are 



a \ _ f 1/24 (-2 + 20a + 5(3) 
P J ~ \ 1/6 (2 + 4a + /?) 



(70) 



and correspond to the states pa,f3- The quantities needed for calculating C are 

Pa,p-P7}s = -;^(4a + 2-5/3)(t/i-f/2) , (71) 

\\Pa,f^ - pT}^\\ = ^ (-4a- 2 + 5/3) , (72) 

{Pa,^,PT}^-Pa,^) = ^(4a + 2-5/?), (73) 
so that operator C is expressed by 

C = ^ (21 + f/i - f/2) . (74) 



The check of the first of conditions (1121) for an entanglement witness gives, unsurprisingly, 

(p-^,^) = -1= (4a + 2-5/3) < 0, (75) 

since 4a < 5/3 — 2 , Eq. flB^ . For the second test we use the fact that operator C ( 1711) is of 
the form (162!) and thus, according to Lemma [3l we obtain 

(p, C*) > \/peS. (76) 

Therefore, C ([71]) is indeed an entanglement witness and the states pa,i3 are the nearest 
separable ones Pa,i3 = Po; a,/3 to the entangled two-parameter states (l50l) of Region II. 
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Finally, for the HS measure of these states, we find 

T/ ent N _ II _,cnt || ^ 

\Pa,l3) — \\P0;a,l3 - Pa,l3\\ 



6^2 



[-Aa -2 + 5/3) 



(77) 



Another way to arrive at the nearest separable states for the two-parameter states is to 



calculate the nearest PPT states with the method of Ref. j28j first and then check if the 
gained states are separable. If we do so we obtain for the nearest PPT states the states 
P/3 and Pa,/3 we found with our "guess method", we know from Ref. 
separable and therefore they have to be the nearest separable states. 



15 these states are 



D. Three parameter entangled states and bound entanglement — qutrits 

1. Detecting bound entanglement with geometric entanglement witnesses 

As we already mentioned, the PPT-criterion 0, |5| is a necessary criterion for separability 
(and sufficient for 2 x 2 or 2 x 3 dimensional Hilbert spaces). A separable state has to 
stay positive semidefinite under partial transposition. Thus, if a density matrix becomes 
indefinite under partial transposition, i.e. one or more eigenvalues are negative, it has to be 
entangled and we call it a NPT entangled state. But there exist entangled states that remain 
positive semidefinite - PPT entangled states - these are called bound entangled states, since 
they cannot be distilled to a maximally entangled state [t', ^ . 

Let us consider states on a di x d2 dimensional Hilbert space, D := did2- The set of 
all PPT states P is convex and compact and contains the set of separable states. Thus 
in Eq. f|TT]) the nearest separable state po can be replaced by the nearest PPT state tq for 
which the minimal distance to the set of PPT states is attained, 

min||r — pII = ||ro — p||. (78) 

If p is a NPT entangled state pnpt and tq the nearest PPT state, then the operator 

AppT := To - Pnpt - (t"o, - pnpt) Id (79) 

defines a tangent hyperplane to the set P for the same geometric reasons as operator f[T^ 
and has to be an entanglement witness since P D S (for convenience we do not normalize 
( I79|) since it does not affect the followin g ca lculations) . The nearest PPT state tq can be 
found using the method provided in Ref. |28|]. In principle, the entanglement of pnpt can be 
measured in experiments that should verify TrAppx Pnpt < 0. If the state tq is separable, it 
has to be the nearest separable state po since the operator (!79|l defines a tangent hyperplane 
to the set of separable states. Therefore, in this case, Appx is an optimal entanglement 
witness, Appx = ^opt, and the HS measure of entanglement can be readily obtained. If tq 
is not separable, that is PPT and entangled, it has to be a bound entangled state. 

Unfortunately, it is not trivial to check if the state Tq is separable or not. As already 
mentioned, it is hard to find evidence of separability, but it might be easier to reveal bound 
entanglement, not only for the state tq but for a whole family of states. A method to detect 
bound entangled states we are going to present now. 

Consider any PPT state pppT and the family of states px that lie on the line between 
Pppt and the maximally mixed (and of course separable) state ;^1d, 

Pa := A PPPT + Id- (80) 
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FIG. 4: Sketch of the presented method to detect bound entanglement with the geometric entangle- 
ment witness C\. The dashed line indicates the detected bound entangled states p\ (Amin < A < 1), 
states p\ with < A < Amin can be separable or bound entangled (straight line). 



We can construct an operator C\ in the following way: 

C\ = p\ - pppt - (pa,Pa - Pppt) Id • (81) 

If we can show that for some Amin < 1 we have Ti pC\^^^ > for all p E S, C\^^^ is an 
entanglement witness (due to the construction of C\ the condition Tr p\Cx < is already 
satisfied) and therefore pppT and all states p\ with Amin < A < 1 are bound entangled (see 



Fig. Hj). In Ref. [30|] a similar approach is used to identify bound entangled states in the 



context of the robustness of entanglement. 

2. Application of the method to detect hound entanglement 

Let us now introduce the following family of three-parameter two-qutrit states: 

1 — (y. — B — ^ 3 1 

Pa,/3,7 := ^ 1 + «^oo + ^ {Pio + P20) + ^ (i^oi + Pu + P21) , (82) 

where the parameters are constrained by the positivity requirement Pa,/3,7 > 0, 

«<^/3 + l-7, «<-/5 + l-7, 

a<-/3 + l + 27, a>f-^ + ^7. (83) 



States (l82l) lie again in the magic simplex and for 7 = we come back to the states (150|) 
considered before. However, for 7 7^ it is not trivial to find the nearest separable states 
since the PPT states do not necessarily coincide with the separable states. But we can use 



our geometric entanglement witness (|HT1) to detect bound entanglement (see Refs. |3ll. |32|). 

We start with the following one-parameter family of two-qutrit states that was introduced 
in Ref. Q: 

Pfc = ^|0'+>(0^| + + 0<6<5, (84) 
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FIG. 5: The states Pa,p,'y represent a pyramid with a triangular base due to the positivity constraints 
(j83p and the PPT points with constraints (j89p form a cone which overlaps with the pyramid. The 
bound entangled and separable states lie in the intersection region. The black dot indicates the 
maximally mixed state and the origin of the non-orthogonal coordinate axes. 



where 



a+ := i (|01) (Oil + |12) (12| + |20) (20| ) , (85) 
a := i (|10) (10| + |21) (21| + |02) (02| ) . (86) 



Let us call this family of states Horodecki states. Interestingly, the states fl8^ are part of 
the three-parameter family fl82|) . namely 



6 - 6 ^ 2b 5-26 
Pb = Pa,f3,^ with a = (3 = - — , 7 = — — , (87) 

and thus lie in the magic simplex. Testing the partial transposition, we find that the 
Horodecki states (18^ split the states in the following way: for < 6 < 1 they are NPT, for 
1 < b < 4 PPT and for 4 < 6 < 5 NPT. In Ref. [9|, it is shown that the states are separable 
for 2 < 6 < 3 and bound entangled for 3 < 6 < 4 . In our case, it is more convenient to use 
7 as the parameter of the Horodecki states. Using Eq. (IHTl) we express b in terms of 7 and 
obtain 

- 1+7 ^ -5 + 77 

Pb = Pa,p,'y With a = —^,(J = — — — ,7. (88) 

The geometry of the three-parameter family of states Pa,i3,-y as part of the magic simplex 
we show in Fig. [5l in particular, the states Pa,i3,-y with positivity requirement (l83l) and the 
PPT states which are constrained by 

^ 1 1 



a 



< ^ (-2 + 11/3 + 3Va) , « < ^ (-2 + 11/? - 3yA) , (89) 
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where A = 4 + 9/3^ + 47 — 77^ — 6/9(2 + 7). The states Pa,/?,^ form due to the positivity 
constraints (!83l) a pyramid with triangular base and the PPT points due to the constraints 
( !89|) a cone. Both objects are quite symmetric and overlap with each other in a way shown 
in Fig. O In the intersection region lie the bound entangled and separable states. 

Application of the method. We now want to apply the method to detect bound entanglement 
of the three-parameter family (l82l) . The idea is to choose PPT starting points on the 
boundary plane a = |/9 + l— 7of the positivity pyramid, on the Horodecki line and in 
a region close to this line, and shift the operators Cx along the parameterized lines that 
connect the starting points with the maximally mixed states. If we can show that Cx is an 
entanglement witness until a certain Xmin, all states px flHOl) with 1 < A < A^m are PPT 
entangled. 

We parameterize our "starting states" on the boundary plane by 

Pplane = Pa,/3,7 With I a = , (j = — , 7 1 , e G R , (90) 

where we introduced an additional parameter e to account for the deviation from the line 
within the boundary plane. 

Depending on 7 and e the operator C^,e,A ( IHTj) has the following form: 



C^,x = Pa - Ppi - {px,px-ppi) 1 = a(21 + ciUi + C2U2 + clU^'' 
d 



with a = -^A(l-A), d = 1 + 87^ + 3e(2 + e)/7 , 
36 

_ 4(2 + e) _ 2(l-7v/37.-3e) 

- '^wT' - ux • ^^^^ 

The operators f/i, f/g, f/g-^ are defined by Eq. fl52l) and the family of states px by 

1 - A 

Pa = A Pplane H — 1 • (92) 

y 

We want to find the minimal A, denoted by Amin, depending on the parameters 7 and e , 
such that all states on the line (l92l) are bound entangled for Amin < A < 1. To accomplish 
this, we define the functions 

5(1(7, e. A) := |ci| and 5(2(7,6, A) := jca] , (93) 

then Amin is attained at max{(7i(7, e. A), 5(2(7, A)} = 1 (recall Lemma[3]). Bound entangle- 
ment can be found in a region where Amin < 1 and the starting points of the lines fl92|) are 
PPT states. That means, e and 7 are chosen such that the starting points are PPT states and 
the corresponding line allows a Amin < 1- The parameter e is bounded by —1/4 < e < 1/3, 
where the lower bound is reached for Amin — > 1 at I7I = 1/4 and the upper bound arises from 
the boundary of the PPT states at 7 = 0. For every e in this interval, we have an interval 
of I7I where bound entangled states are located. 

More precisely, in the interval for the e parameter —1/4 < e < eo with eo = 
(8 - 7(2/(-5 + y29))i/3 ^ 7(2/(-5 + v^))-^/^^ 3 ^ _o.o3 the parameter 7 is con- 
fined by Vl-2e + 3e2/v^ < I7I < \/7-Qe- 3e2 - 2 (1 - 48e - 12e2)V2 /^, under 
the constraint Amin < 1- For the remaining interval eo < e < 1/3, we get the bounds 
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FIG. 6: Plot of the parameter range 7 versus e of starting points of the hnes ([92]) such that bound 
entangled states can be detected. 

< I7I < \/9 — 26e — 3e2 / 7 , where the lower bound is again constrained 
by Amin < 1 and the upper one by the PPT condition. A plot of the allowed values of e and 7 
for the starting points on the boundary plane is depicted in Fig. [61 We have equality of the 
coefficient functions (71 = (72 for I7I = a/15 + 22e — Se^ / 7a/3 =: 70 , (/i > (72 for I7I < 70 and 
S'l < ^72 for I7I > 70) where I7I is always restricted to the allowed range described above. As 
mentioned, Amin is gained from the condition max{5fi(7, e, A), (72(7, e, A)} = 1 for particular 
values of 7 and e. The total minimum Aj^*^ is finally reached at 



which is significantly below the value 1 so that the resulting volume of bound entangled states 
is remarkably large. The total minimum ( [Ml) is attained at e = (TvTS — 25) / 2 ~ 0.12 and 
I7I ~ 0.35 . The whole line of states px i M2\i within the interval X^^l^ < A < 1 is found to 
be bound entangled. The volume of the detected bound entangled states we have visualized 
in Fig. [3 For 7 = no bound entanglement occurs, as discussed in Sec. IIII C[ which is 
represented in Fig. [71 at the meeting point of the two bound entangled regions. 



In the last section, we gave a strict application of our method to detect bound entangled 
states, where we recognized bound entangled states on the parameterized lines p\ ( 180]) only. 
The involved entanglement witnesses, however, are able to detect the entanglement of all 
states on one side of the corresponding plane, not only on the lines. Therefore even larger 
regions of bound entanglement can be identified for the three-parameter family ([82]) . which 
is described in detail in Ref. j32| . 

In this subsection, we want to show that for the three-parameter family ([52]) our method 
detects the same bound entangled states as the realignment criterion does, and furthermore 
allows for a construction of the shape of the separable states of the three-parameter family, 
so that we are able to fully determine the entanglement properties of this family of states. 

The realignment criterion is a necessary criterion of separability and says that for any 
separable state the sum of the singular values Si of a realigned density matrix has to be 




(94) 



3. More bound entangled states and the shape of separable states 
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FIG. 7: Regions of detected bound entangled states within the pyramid represented by the three- 
parameter states Pa,i3,'y (|82|) . The dot represents the maximally mixed state, the Horodecki states 
are represented by the line through the boundary plane from which the regions of bound entan- 
glement emerge. In the magnified picture on the right hand side the viewpoint is altered a bit to 
gain a better visibility. 



smaller than or equal to one, 



TrWorjjfjR < 1 



(95) 



where {pij^ki)^ '■= Pikji (for details see Refs. [33|, ISJ, l35|, l36|]). States that violate the criterion 
have to be entangled, states that satisfy it can be entangled or separable. 
In our case of the three-parameter family fl82l) we obtain the constraints 
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(6 + 11/5-7- Ai) 
(6 + ll/?-7 + Ai) 
(-6 + 11/? - 7 - A2) 
(-6 + 11/3 - 7 + A2) 



from the realignment criterion, where 



Ai := V4 + 36/5 + 8I/52 - 127 - 54/37 + 3372 and 
A2 := ^/4: - 36/5 + 81/52 + I27 - 54/^7 + 337^. 



(96) 
(97) 
(98) 
(99) 



(100) 



Only constraint (!96|) is violated by some PPT states, which thus have to be bound entangled. 
The PPT entangled states exposed by the realignment criterion are therefore concentrated 
in the region confined by the constraints 



a < -/? + 1 - 7, a 



< ^(-2 + ll/5 + 3v^), a> ^(6 + 11/5-7- Ai) 



101) 
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All PPT entangled states of Eq. fllOip can also be detected using Lemma [31 To see this, we 
construct tangent planes onto the surface of the function 

a = 4(6 + ll/?-7-Ai) (102) 

from the realignment criterion ([HH]), where we use orthogonal coordinates. In this way, we 
can assign geometric operators to the tangential planes by choosing points a inside the 
planes and points h outside the planes such that a — 6 is orthogonal to the planes. Since the 
Euclidean geometry of our picture is isomorph to the Hilbert- Schmidt geometry, the points 
a and h correspond to states pa and pf, and we can construct an operator accordingly, 

Cre = Pa-Pb - {pa, Pa - Pb)^9 ■ (103) 

Decomposed into the Weyl operator basis, the operators fll03p that correspond to tangent 
planes in points (at, A^Tt) ^ where at is a function of (3t and 7t, given by the realignment 
function (11021) - are 

= a (2 1 - Ui + cUi + c*U^^) , with 

a = 4 (-2 - 9A + 37t + 3A,) , 
3d 

_ H + (-2 - 9A + 37t) Ac + y37t (2 + 9A - 37t + 3Ae) t 
(2 + 9A)^-6(2 + 9A)7t + 36 7t^ 

A, := ^4 + 36 + - 127t - 54A7t + 337? ■ (104) 

The absolute values of the coefficients c and c* in Eq. (11041) are 1, |c| = |c*| = 1, and 
therefore, according to Lemma [H the operators Cre are entanglement witnesses that detect 
the entanglement of all states "above" the corresponding planes, thus also the bound 
entangled states in the region of Eq. (llOip . The detected bound entangled region is depicted 
in Fig. [HI Naturally, now the question arises if all the states Pa,f3,-y (j82l) that satisfy both the 
PPT and the realignment criterion are separable. This cannot be seen using the two criteria 
alone, since they are not sufficient criteria for separability. But we can apply a method of 
shifting operators along parameterized lines in the other direction: by constructing a kernel 
polytope of necessarily separable states and assigning operators to the boundary planes 
of this polytope, we can shift those operators outside until they become entanglement 
witnesses, which can be shown numerically (unfortunately Lemma [3] does not help here). In 
this way, one can step by step reconstruct the shape of the separable states, which indeed 
is the set given by the states that satisfy the PPT and the realignment criterion (see dark 
shape of the two intersecting cones in Fig. [H]); for details see Ref. [13]. Hence, we see that 
the presented method of "shifting" operators along parameterized lines does not only help 
to detect bound entanglement, but also to construct the shape of the separable states. 



IV. CONCLUSION 

We discuss the geometric aspects of entanglement for density matrices within a simplex 
formed by Bell states. We use entanglement witnesses in order to quantify entanglement 
and detect in case of qutrits bound entangled states in specific instances. 
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FIG. 8: Illustration of the PPT and realignment criteria. The familiar PPT cone (here pictured 
only inside the positivity pyramid) is intersected by a cone formed by the realignment constraint 
(j96p and reveals a region of bound entangled states (translucent yellow region) that can also be 
detected by Lemma El 

We demonstrate the geometry of separability and entanglement in case of qubits by 
choosing so-called two-parameter states, Eq. (fT5i) . i.e., planes in the simplex, a tetrahedron 
(see Fig. [2]). These states reflect already the underlying geometry of the Hilbert space and 
they are chosen with regard to the description of qutrit states, a generalization into higher 
dimensions. To a given entangled state we determine the nearest separable state, calculate 
the corresponding entanglement witness and the Hilbert-Schmidt measure in the relevant 
Regions I and II (see Fig. [1]). 

In case of qutrits it is quite illustrative to demonstrate the geometry of separability and 
entanglement in terms of two-parameter states (150|) . These states represent a plane in the 
eight-dimensional simplex formed by the nine Bell states, the magic simplex, and are easy 
to construct within the Weyl operator basis. Due to the high symmetry of this simplex we 
may restrict ourselves to a certain mixture of Bell states, Eq. (130]) . which lie on a line in a 
phase space description. This line exhibits the same geometry as other lines. Within the 
Weyl operator basis it is quite easy to find the Bloch vector form fl5T|) of the two-parameter 
states. It enables us to find in regions I and II (see Fig. [3]) the nearest separable state to 
a given entangled state and the corresponding entanglement witness. The easy calculation 
of the Hilbert-Schmidt measure of entanglement is a great advantage and its result of high 
interest. Other entanglement measures, like the entanglement of formation, are much harder 
to be calculated in this higher dimensional case. 

We present a method to find analytically bound entangled states by using entanglement 
witnesses. These witnesses are constructed geometrically, Eq. flSTl) . in quite the same way 
as for the detection of the nearest separable state, Eq. f[T^ . We show that the Horodecki 
states (18^ can be described by the family of three-parameter states (l82l) and are therefore 
part of the magic simplex. Geometrically, they form a line on the boundary of the pyramid 
represented by the three-parameter states (!82|) (see Figs. [7|and[8]). We apply our method 
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to find regions of bound entangled states within the pyramid of states (l82l) (see Fig. [7j). 
Even when restricting ourselves strictly to consider the detected bound entangled states on 
the parameterized lines (IHOj) only, we find large regions of bound entanglement. Employing 
the realignment criterion of separability, we can reveal larger regions of bound entanglement 
that are also detected by Lemma [3l Finally, we can apply our method of shifting operators 
along parameterized lines together with numerical calculations to show that there do not 
exist more PPT entangled states for the three-parameter family. Hence, the shape of the 
separable states can be constructed. 

When decomposing density matrices into operator bases the Weyl operator basis is the 
optimal one for all our calculations. The reason is that entangled states - the maximally 
entangled Bell states - are in fact easily constructed by unitary operators a la Weyl, see 
Eq. die]). 
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